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Abstract 

The non-Abelian Kubo formula is derived from the kinetic theory. That 
expression is compared with the one obtained usingthe eiknoal for a Chern- 
Simons theory. The multiple time-scale method is used to study the non- 
Abelian Kubo formula, and the damping rate for longitudinal color waves is 
computed. 
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1. Introduction 

Substantive progress is made in the study of non-equilibrium quark-gluon plasma in recent 
years Q. The Kubo formula is a fundamental tool for describing a non-equilibrium system. 
Hence it is important to get the Kubo formula of a non- Abelian system. In Ref. [2] , meaningful 
results are given. In this paper, we derive a form of the non- Abelian generation of the Kubo 
formula from kinetic theory directly, and the multiple time-scale method is used to investigate 
it. To show the intention of our work concretely, we first give the result in Ref. [2]. Using 
the relation between the current and the generating fuctional, i.e., = —8T[A]/ 8A^(x), 

and the Yang-Mills equation, a form of non- Abelian Kubo formula is obtained with the color 
current being expressed in terms of the eikonal for a chern-Simons theory ||: 

[D v F^(x)] a =j%(x), (1) 
where D u = d u - ig[A u (x), • • •], 

fa(x) = ^ / d 4 ke~^^(k) = J d*ke-* k -*jr^ a (k), (2) 

and jfi a (k) is expressed in terms of powers of A^, 



n=l 



fciH Vk n =k 



\x{{I a v»)(A{k{) ■ v)}8 nl 



g n -Hv\ (IaV^Aih) •«)... (A(k n ) • v) }F(h, ■ ■■ , k n ) 



(3) 



F(ki, ■ ■ ■ ,k n ) = 7= ^TF =1 F Ql - \r- ~\ T =T> ( 4 ) 

^ (9o - 9i)(9i - ft) ■ ■ ■ (9i-i ~ ft)(9i+i ~ 9i) ■ ■■ (9n - ft) 

i 

where v — (1, v) is a four light-like vector, and v' = (1, — v), q~o = 0, ft = £ ' v + 

3=0 

i 

q L = £ {kj ■ v'). Nf being the number of flavors for massless quarks, I a is the SU(iV c ) group 
j=0 

generator in the fundamental representation. d£l denotes integration over angular direction 
of the unit vector v. In the above expression, the relation J dflf[A ■ v] = j dSlf[A ■ v'] has 
been taken consider of, and we have defined, 



The non-Abelian Kubo formula describes Landau damping of fields in a quark-gluon 
plasma as well as Debye screening and progration of plasma waves in a gauge covariant way. 
Finding solutions to the non-Abelian Kubo formula is appealing, since such solutions give 
collective effects of the quark-gluon plasma at high temperature. R. Jackiw, Q.Liu and C. 
Lucchesi investigate the static non-Abelian Kubo formula and show that there are no hard 
thermal solitonsQ]. However, the above expression of the Kubo formula is quite complicated, 
this makes it much more difficult to use it to analysis relativistic physical problems. Recently, 
P. F. Kelly, Q.Liu C. Lucchesi and C. Manuel have shown that the generating functional of 
hard thermal loops can be derived from classical transport theory]^], and it has been shown 
that the solution of the kinetic equations taking into account the non-Abelian and nonlinear 
effects can be obtained using the multiple time-scale method ||. The multiple time-scale 
method is one of the most powerful asymptotic methods in solving nonlinear equations, and 
it is proved successful in studying nonlinear properties of electromagnetic plasmas [7]. This 
gives us a chance of investigating the non-Abelian Kubo formula from kinetic theory. In 
this paper, the non-Abelian Kubo formula will be derived from the kinetic theory, then the 
multiple time-scale method will be used to study it and damping rate will be calculated. 

This paper is organized as follows: In section 2, the non-Abelian Kubo formula is derived 
from the kinetic theory and the identity of it with the form obtained using the eiknoal for a 
Chern-Simons theory is shown. In section 3, the multiple time-scale method is introduced to 
study the non-Abelian generation of the Kubo formula. In section 4, using the multiple time- 
scale method, the damping rate due to non-Abelian and nonlinear nature of QGP cigenwaves 
is calculated in the temporal gauge for longitudinal color waves, and the result is compared 
with the hard thermal result in thermal QCD. 
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2. Non-Abelian Kubo formula derived from the kinetic 
approach 

In this section, we will derive the non-Abclian Kubo formula from the kinetic approach. 
The phase space of the distribution function of colored particles is (x,p, Q), where Q a , a — 
1, 2 • • ■ TV 2 — 1 are color charges, and the use of the color charge as phase-space coordinates 
is justified recently [5]. Taking cosider of physical constraints on the phase space volume 
elemeant dxdPdQ, the momentum and the color charges measure 

dP=-^m P °)5(p 2 - m i), (6) 

dQ = d 8 QS(Q a Q a - q 2 )8(d abc Q a Q b Q c - q 3 ), (7) 

where the constants qi and q 3 fix the values of the Casimir invariants and d a bc are the totally 
symmetric group constants. 

When the classical transport equations for colored particles is expanded consistently in 
the coupling constant g, at the leading order, the transport equations for the distribution 
functon f(x,p,Q) is reduced to(neglecting spin effects) [5] 

P"fr. ~9f abc AlQ c ]f(x lP ,Q) = - gQa F^(x)-^-f^(p), (8) 

where distribution functions f^°\p) in the above equation is the equilibrium distribution 
function in the absence of a net color color field, 

fW(p)==Cn B ,F, (9) 

where C is a normalization constant, ub,f are Fermi-Dirac and Bose-Einstein equilibrium 
distribution function respectively 

n F = (e^ u + n B = {e pp U (10) 

where /3 is the inverse of the system temperature T, is the local four velocity. 
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The Yang-Mills mean field equation coupled with the transport equation is 
[D v F^(x)] a =j£(x), (11) 
where color current j u generated by quarks and gluons is expressed as, 

j"W= E E / dPdQQ a f(x,p,Q) (12) 

species hclicites 

Using the above equation, the constraint on the first order current demanding by the leading 
order transport equation can be obtained^] 

[p • DJ»(x,p)] a = 2g 2 p^F? (x)-^[N f fW(p) + JV c G<°>(p)], (13) 

apo 

where the field tensor is only related to F v q. 

The constraint equation (13) can lead the generating functional of hard thermal loops ||. 
This means that Eq.(13) determines the color current generated by hard thermal loops and 
can be used to study the response of quark-gluon plasma to fields consistently. In the 
following discussion, we will use the constraint equation (13) to get the expression for the 
color current. 

In the momentum space the constraint equation ( [l3]) yields, 
(p-k)J^k,p)+tgr bc £ (p.A(k%J?(k",p) 

k'+k"=k 

= 2flV [(p • fc)ffo wp v ]Al{k)-^[N c G^{p) + JV)/W(p)] 

+2i fl Vy E f abc Al(ki)A c (k 2 )^[N c G^(p)+N f f^(p)}. (14) 

k 1 +k 2 =k 

The color current can be expressed as a series of powers in A^, and in the above equation 
there is no constraint on the field potential A^, thus Eq.(|l4|) can be reduced to the equations 
of every power in A^ . After integrating over the momentum space one can get 

2 fdn CJV^V" 



#.(*) = -mlKti - / ^ v ; +iQ+ ] (15) 



(v ■k + iO+)(v ■ ki +i0+) 

ki+k 2 =k 

v»v x v v A b v (k 2 )Al(h), (16) 



jH <a (k) = -igf abc [dSl 1 J2 {v-A{k')) b] ^\ c (k") n>2, (17) 

J k' + k" = k 

where m 2 D = ±(N C + \N f )g 2 T 2 , j M n _ liQ (A:) = / dnj^\ a (k) and to avoid the poles in the 
above integrand, the retard condition, i.e., replacing uj by lo + i0 + , has been imposed. 

The expression for the current to any power of can be obtained easily using the 
above equation. Further, j£ can be expressed as 

+(-ig) n - 1 m 2 D f — r a " fc "- 1 / b - ia "- lb "- 2 ---/ b2Q2bl 

^ n k 1 + -+k n =k 

(vk + iO+)(v ■ (fc - k n ) + i0+) ■■■(v(k-k n k 2 ) + i0+) 

(A(k n ) ■ v) an ■ ■ ■ (A(k 2 ) ■ v) a2 {A{k!) ■ v) bl . (18) 

Using Eq.(2) and the above equation to express the color current in the Yang-Mills mean 
field equation (11), the non-Abelian Kubo formula is derived from the kinetic theory. The 
above expression for the color current is relatively simple, plain, and the non-Abelian Kubo 
formula expressed by it is easier for being studied. In the case when the mean field is weak, 
the non-Abelian Kubo formula can be solved by the multiple time-scale method. 

Now we turn to show that the expression for the color current obtained from the kinetic 
theory is equivalent to the one obtained using the eiknoal for a Chern-Simons theory. 

In the kinetic theory, the expression for the color current is determined by Eqs.( |l^- |l7j ). 
One can check that j'f a (k) and j% a (k) in Eq.( |l5|) and Eq.( |l6|) are equivalent to those got 
from Eqs. (3- ^) . Thus to show the identity of these two expressions, one only need to show 
j% a in Eqs. (3- ^) satisfy Eq.( [l?]). Now we turn to show this. 

Using the relation 

1 1 1 



qi(q n - qi) q%q n q n {q n - qi) ' 

And taking note of the delta function in the definition (5), one can get from Eq.(4), 



(19) 



F(h, ■ ■ ■ , kn) = - ]_ . n+ [F(fcl, • • • , kn-!) - F(k 2 , • • • , k n )}. (20) 

v ■ k + ilr 



Then using the invariance of j£ a when fei, • • k n are exchanged, we obtained from Eq.(3) 
1 7V T 1 ^ f f 

(A(fcx) •«)... (A(fc n _i) • (fci, • • • , k n -i){A(k n ) ■ v) b 
= l d ^ v . k \ l0+ E {vA{k')) b f n \ c {k"). (21) 

J k'+k"=k 

This concludes our proof of the identity of the expression of the non-Abelian Kubo formula 
from kinetic theory and the one obtained using the eiknoal for a Chern-Simons theory. 

In the following section, we will investigate the non-Abelian Kubo formula obtained 
from the kinetic theory, i.e., Eq.(ll) coupled with Eq.(2) and Eq.( |l8| ), using the multiple 
time-scale method. 



3. The multiple time-scale method 

The multiple time-scale method is a kind of asymptotic methods of solving nonlinear equa- 
tions. In plasma theory, it is used to analysis the nonlinear effects of a electromagnetic 
plasma[7]. Generally speaking, when the field strength is not very strong, the mean field 
equation may be solved by perturbation methods. However, As is known, the nonlinear 
wave equation in a plasma can not be solved by naive perturbation analysis, which may 
yield singularity in the solutions by virtue of the resonance in the relevant equations. This 
difficulty may be overcome if it is solved by means of the multiple time-scale perturbation 
approach[7,8]. In this method, first of all, the independent time variable t is extended to 
many time variables, each with different scale, by introducing the new time variables [7] 

t (n) = K n t, n = 0,l,2,---,JV, (22) 

where k is a dimensionless parameter introduced to denote the order of a small quantity. 
Then, every function dependent on t, such as Y[t], is expanded into an asymptotic series of 
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the form 

N 



Y[t, K } = J2^Y (n) [t ( "\t il \---,t (N) }+0(K N+1 ). (23) 

n=0 

Every function dependent on t is assumed to be a function of the multiple time, so the 
time derivative should be replaced by 

Usually, one prefers to discuss the nonlinear wave equations in momentum space. The 
multiple time-scale procedure can also be used there by introducing the multiple scale Fourier 
transformation 

poo />oo 



0J K ' , ■ ■ ■ , W v 

N \ N 



/oo poo 
.../ r(,<V . . 
-co J —oo 

(N v iv 

-i^wWtW JJdwW. (25) 
n=0 ' n=0 

Using this transformation, the multiple time-scale method can be used in momentum space 
simply by taking the following expansions instead of Eqs.(22-24)[8], 

N 

w = ^A<"', (26) 

n=0 

N 

Y[u, k] = Y, K " y(n) ^ (0) > w(1) > ■ • • ' w(Ar) ] + °^ N+1 )- (27) 

In this paper, we will use the multiple time-scale method to study the non-Abelian Kubo 
formula in momentum space. The corresponding equation of Eq.(ll) in momentum space is 
obtained by using 

A^x) = J -0^ e - ik -*A»(k). (28) 
The non-Abelian Kubo formula in the momentum space is 

- (k 2 g"»-k»k»)A v (k)-g ]T {k v + k%)[A l/ (k 1 ),A»(k 2 )] 

ki+k 2 = k 

+ g ]T {K)[A v {k 1 \A v {h 2 )]+g i E [^(k^iAMA^ki)}} 

k\-\-k2—k k±-\-k2—k fc 3 +fe 4 =fe2 

= f{k), (29) 



where j^(k) is expressed as -Eg. (18). 

As we have discussed in the beginning of this section, when the mean field is not very 
strong, the non-Abelian Kubo formula can be solved by the multiple time-scale method. If 
we work in the momentum space, following Eq.( p6| ) and Eq.( |2^ ), the frequency and the 
color field potential are expanded into asymptotic series of forms as 

= K (JV iW W, w (2) 5 w (3) t . . .) + K 2 A {2) (w (0) ^ w (l) ; w (2) ^ w (3) ^ . . .) 

+ k* A® (J® , u^, u<®, •••) + •••■ (31) 

Inserting the expansions of all the variables into the nonlinear equation, and equating the 
coefficients of corresponding powers of k in the two sides of this equation, the perturbation 
expansions of the nonlinear equation at every order are given. The nonlinear equation then 
can be discussed iteratively. 

Now, we first discuss the physical meaning of the above expansions given by the multiple 
time-scale method. In the linear response approximation, any perturbation in a plasma 
may be taken to be a superposition of independent eigenwaves whose dispersion relation is 
determined by the solution to the non-Abelian Kubo formula in the linear approximation. 
When the nonlinear terms in the non-Abelian Kubo formula are retained, the nonlinear 
interactions of the eigenwaves will be taken into account. The nonlinear interactions of the 
eigenwaves will change the dispersion relation obtained in the linear approximation, i.e., 
frequency corrections expressed by Eq.( |30| ) will be given to the eigenfrequency 

cj(°) correspondingly. If the imaginary parts of the frequency corrections for eigenmodes 
are not vanishing, there will be instablity or damping for the eigenmodes. In next section, 
we will calculate the damping rate due to interactions of QGP eigenwaves in the temporal 
gauge. 
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4. Nonlinear damping rate in the temporal gauge 

The non-Abelian Kubo formula determines the dispersion relation of the eigenwaves in 
QGP. The damping rate of plasma eigenwaves is determined by the imaginary part of the 
dispersion relation. In this section we will calculate the damping rate for longitudinal color 
eigenwaves in the long-wavelength limit from the non-Abelian Kubo formula. 

In this section, for convenience we will work in the temporal gauge, A® a = 0. The relation 
between the color electric field E and the vector potential in that gauge is simple, 

E\{x) = -dA\{x)/dt. (32) 

The non-Abelian Kubo formula in the momentum space and in the temporal gauge is 

-(w 2 -K 2 )A h (k)~k l k h A\k)+g k i [A i (k 1 ),A h (k 2 )} 

k% + k 2 =k 

+g ]T *&[A i (k 1 ),A h (k a )]-g £ A* fa)] 

k\-\-k2 — k k\ J \-k2 — k 

+9 2 E E [A i (k 1 ),[A i (k s ),A h (k A )]]=j h (k), (33) 

where K = |k| and j h (k) is expressed as Eq.(18). 

Though the method used in this paper is general, for simplicity of calculation, we suppose 
that there is only longitudinal vector potential, i.e., the vector potential is parallel to the 
wave vector, 

A\k) = k l A{k)/K. (34) 

Inserting the expansions ( |3l]) into Eq.(33) and Eq.(18), for longitudinal field, the 
first order non-Abelian Kubo formula is, 

- M°) 2 - K*)A^{k) - *«*MW(Jb) = -ml J (35) 

The above non-Abelian Kubo formula can be reduced to 

e(wW,k)i4W(A) =0, (36) 
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where e is the color electric permeability, 

^ (0) >k) = i + |f [i - ^Nf^l ~ W ~ - (0) ))]- (37) 

where oj 2 = (2iV e + N s )g 2 T 2 /18. 

From the above equation, the solution to the first order non-Abelian Kubo formula can 
be got 

= -^^o[^<5M 0) - w£) + e **E«y( w C0) + (38) 

where -E£ and <^>£ are the initial amplitude and phase of the oscillation respectively (a is 
introduced to label different eigenwaves). The frequencies and the wave vectors must satisfy 
the dispersion relation, 

£(J°),k) = 0. (39) 

Eq. ( |3^ ) , which is quite similar to the dispersion relation for the abelian-like plasma eigen- 
waves, is the first order dispersion relation for the color eigenwaves in a QGP. It agrees with 
the leading order result using the finite temperature QCD||. It has been shown by U. Heinz 
that as the eigenwaves in QGP are always time-like, lu^ 2 > K 2 , the color eigenwaves in a 
QGP do not undergo Landau damping in the linear approximation[10]. 

Using the results of the first order equation, the second order equation now can be 
considered. The second order non-Abelian Kubo formula for the longitudinal field is 

-^\(^\\,)A( 2 \k)^-2u^^A^(k)^ 

k±+k 2 =k 



' D 



^ (0 W>( fc) - vk 



4tt K a w v ■ k + i0+ 



-2gm 2 D g J 



dn v ki v k 2 v h uf ] 

4^ KxKi (v ■ k + iO+)(v ■ h + i0+) 



trj^^ 1 )^),^ 1 )^)]}, (40) 

where the right side of the above equation comes from the second order contribution of j£ 1 
and j% 2 in Eqs. (15-16). 
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Eq.( 40) discribes the nonlinear effects owing to the three- wave processes which are 



the merging of two eigenwaves with eigenfrequencies 0J^ a , * into a third wave, or the 
corresponding inverse processes. The appearance of the term [A^(ki), (A^)] in Eq.( p0| ) 
indicates that the non-Abelian nature of the mean color field has been taken into account 
and the interaction of the QGP eigenwaves. 

Now, we consider the nonlinear effects due to the three- wave interactions. We will see 
that the nonlinear effects owing to the three-wave interactions on the eigenfrequency are 
vanishing. In the three-wave processes, when the third wave is also an eigenwave, this kind 
of processes is called three- wave resonance. Taking note of Eq.(5), the following match- 
ing condition representing the conservation of energy should be satisfied for the three-wave 
resonance, 

i4f ±4° 2 )a i=4 o)t , (4i) 

where cr, A, and r in the above equation indicate that the corresponding frequency is the 
eigenfrequency of a certain eigenmode. 

The three- wave resonance processes require that the matching condition Eq. ( ^l]) must 
be satisfied strictly; even a small mismatch of Eq. ( |4TJ) may significantly decrease the effi- 
ciency of the three-wave resonance[ll]. This condition limits the permissible region of the 
three-wave resonance interaction in k space. Indeed for some forms of dispersion relation, the 
three- wave resonance is completely forbidden. It can be proved that for dispersion Eq.( |37| ) 
and Eq.( [39]), where the phase velocity decreases with K [1,10], the three- wave resonance is 
forbidden jl2]. The absence of the three- wave resonance here leads the eigenfrequency cor- 
rection a/ 1 -* to vanish, which is similar to the corresponding situation in an electromagnetic 
plasma [8]. Using this result and taking note of Eq.(38), the second order field potential is 
expressed as 
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-2gm%g x 



ki+k 2 =k 



dn v • k lV • k 2 v ,i uj ( ° ] 



4tt ifiX 2 (v • k + iO+)(v ■ ki + iO+) 
tvilaiA^ih^A^ih)}}. (42) 



From the above discussion, we can see that there is no effect of three-wave resonance 
on the eigenfrequencies. Hence the next order perturbation equation should be considered 
in order to obtain non-vanishing effects of the non-Abelian and nonlinear wave interactions 
on the eigenfrequencies. 

Using the above result Eq.( ^) and a/ 1 - 1 = 0, the third order non-Abelian Kubo formula 
for the longitudinal field is expressed as 

-^°) 2 e(u/°) , k)4 3 ) (fc) - 2^ 2 U°) A™ (fc) 

+ 2 * 2 E E ^^{/.[AW^),^)^),^)^)]]} 

ki-\-k2—kk^-\-k^—k2 



j.. i ;„ j . 7 . i j.. r. . « 



v ■ k v • kj v • k 2 



0+ K K x K 



2 



4 0) ^ (0) , 1 1 /"„,, 1 



«-4 0) +i0+ u-fcj 0) +i0+ J° )2 e(k { 2 ] )J v'-ki 0) +i0+ 

u' • fc£ J + i0+ ^2 ^3 ii4 

o22 4 f Jn 1 v • k v • k 3 v • ki v • k 4 
-2q z m n > ? Q dtt -r ; 

7„ I 7„ I- f. I 7- 7, •* 1 tJ - 1 - ^ 



tr{Z [AW(fc 1 ) J [A«(fc 8 ),A( 1 )(fc 4 )]]}. (43) 



1 , , (0) 

1 LO A 



v ■ kf } +i0+v kf } + 10+ 
We only need to discuss the four-wave resonance processes, as we are interested in the 
effects of wave interactions on the eigenfrequencies. Therefore, k^ " 1 in Eq.( ^ ) satisfies the 
dispersion relation Eq.( ^) and Eq.( ^ ) in the following discussion. In this section, we only 
calculate damping rate being the imaginary part of the eigenfrequency, i.e., 7 = -Imw' 2 '. 



Since only the imaginary part of Eq.( 43 ) is related to the damping rate, it is also sufficient 
to consider the imaginary part of Eq. ( ^ ) 

As the oscillations are developed from random thermal motions, the eigenmodcs can be 
considered completely incoherant to each other. For a color field being in random phase, we 
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have 

(A (1)Q (fc)) = 0. (44) 

where ( } means the average with respect to the random phase of the oscillations. Then, 
using Eq.( |3|), one can obtain 

(A^ a {k)A^ b *(k')) = {2v)H^{k - k')6 b a {A^ 2 )^, (45) 



<^ (1)2 >k. = -^5 [<^ (0) - Uk) + S(cjM + Wk )]/ k; (46) 

UJ\°> 

where Ik characterizes the total intensity of the fluctuation oscillation with frequency Wk 
and — cjk, V is the volume of the plasma. 

As the oscillations are incoherent to each other, the average of product of the field 
potential can be expanded as 

{A^ik^A^i^A^ik^A^ih)) 

= (A^(k 1 )A^(k 2 ))(A^(k 3 )A^(k i )) + (^( 1 )(fci)^( 1 )(fc 3 ))(^ (1) (fc 2 )A( 1 )(fc 4 )) 
+ (AW ( kl ) A W (fc 4 )) (A« (fc 2 )A« (fc 3 )>. (48) 

Multiplying both sides of Eq.(^) by A^ b {k'), then performing the average of the result 
with respect to the random phase using Eq.(||) and Eq.(|f), one can obtain the damping 
rate of the longitudinal eigenwaves in QGP 

j(K) = -Imw (2) 



4tt J (2tt) 4X 1 1 vk(°) +i0+ 

r ( vk 2 vk l 2 1 L>M 

^ K Ki v(kW -k^) + i0+ \.k[ 0) -i0+ vkW+tO+> 



Vk Vvki V-(k-k!) ^°)-^ 0) 

K H K x H |k-ki| )( v . {k (o)_ k (o) ) + l0+ v . k W +lQ+ < 



d£V_ v'-(k-ki) wv'-k v'-ki 1 1 

a— \ m i, i a j? a 7? i 



4tt v |k-ki| K Ki ; e ( w (o) - w ( ),k-ki) (w(0) -^ 0) )2 
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u' • (fc(°) - fci 0) ) + «)+ V • k{ 0) - z0+ «' ' fc (0) + «0+' 



(49) 



Now we perform the integrals in Eq.( ^9| ) in cylindrical coordinates and in plasma 
particle local rest frame. For simplicity of calculation the direction of polar axis is selected to 
be the same as the direction of ki; then the coordinates for k, p, p' arc \t(K, a, (3), p(E p , 9, tp), 
p'(E' , 9', tp') correspondingly. As only the damping rate in the long- wavelength limit is 
calculated numerically later, we will use the approximate relation in the long-wavelength 
region in the following calculation. In the long-wavelength region, the dispersion relation 
Eq.(39) is reduced to the form of [1,10] 

J^=u;1 + \k\ (50) 

For the long-wavelength limit case(-?£T = 0), using the equation p°/(p ■ k' + ip°Q + ) = 
P[1/(lu' — v • k')] — m8(Lo' — v • k'), where v = p/p° and P stands for principal value of the 
function, one can check that the condition for nonvanishing of Eq.(49) is the satisfation of 
the delta function <5[o/°) — lj^ — v ■ (k — ki)]. This delta function restricts the upper bound 
of the integrals over K\ in Eq.(49) to equal gT approximately. We obtain from Eq.( |^) the 
expression of the damping rate for k = modes, 

2™2 r r aT ,.,(0) 



1 d(cos9) p ~d<pS[cos9 - ^ k ^ ] J ™L 

2ir 



f m 2 f 1 f 2lT 
<yCos 2 (j)cos 2 a -\ J d(cos8') J dip 1 



(uj p - Lu[ 0) ) 2 e(uj p - ujf \ -ki) 



(o) v ai C°) 

,,r,Lo\' K 1 cos9 . w - w„ 

cos^coscj)' - 1 + —- L v - 

u v lu^> - K lC os9' lj p - + K lC osB' 

,,(o) ,.(0) >, 

+ ^ ( 51 ) 

where </>, <fr', are the corresponding angles between p and k , p' and k. They can be expressed 
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as 

coscj) = sin6sinacos(cp — 0) + cosOcosa. (52) 

And the expression for cos<^' is similar to the above equation. 

In Eqs.( pl| ), the damping rate is proportional to the wave intensity 7k ■ In an equilibrium 
plasma, the wave intensity is governed by the plasma temperture[7][13], 

/k = 8tt-^—, (53) 
e k w k 

For the long- wavelength case, using Eq.(50), we can obtain 

I k = 4ttT. (55) 

Then, after finishing all the integrals in Eq.( pl| ), the numerical result of the damping 
rate for pure gluon gas is obtained 

7 - 0.21.g 2 T. (56) 

Now, compare the above result with the damping rate for long-wavelength gluons cal- 
culated by the effective perturbation theory based on the resummation of the hard thermal 
loops. The hard thermal result is [14] 

7 - 0.26.g 2 T. (57) 

Note that both of them are obtained approximately, we can conclude that they are coincide. 
This indicates the consistency between the high temperature limit of thermal QCD and the 
non-Abelian Kubo formula derived from the kinetic theory for quarks and gluons. 

In this paper, we start from the non-Abelian Kubo formula which is gauge invariant, 
and we chose a physical gauge to express the gauge invarant equation, then the multiple 
time-scale method is used to expand the equation according the field strength. Our result is 
coincide with the gauge invariant result obtained using the effective perturbation theory. As 
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to the question if the multiple method guarantees the gauge symmtry generally, this need 
further study. 

In conclusion: in this paper we derive the non-Abelian generation of the Kubo formula 
from the kinetic theory and show it is coincide with the form obtained using the eiknoal for 
a Chern-Simons theory. The multiple time-scale method is used to solve the non-Abelian 
Kubo formula obtained in this paper to calculate the damping rate for the eigenwaves in the 
long-wavelength limit numerically 

the result is coincide with that obtained using the effective perturbation theory in the 
finite temperature field theory. This show that the multiple time-scale method is effective 
on solving the non-Abelian generation of the Kubo formula. 
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